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ABSTRACT 

We consider the propagation of cosmic rays in turbulent magnetic fields. We use the models of 
magnetohydrodynamic turbulence that were tested in numerical simulations, in which the turbulence 
is injected on large scale and cascades to small scales. Our attention is focused on the models of 
the strong turbulence, but we also briefly discuss the effects that the weak turbulence and the slab 
Alfvenic perturbations can have. The latter are likely to emerge as a result of instabilities with in 
the cosmic ray fluid itself, e.g., beaming and gyroresonance instabilities of cosmic rays. To describe 
the interaction of cosmic rays with magnetic perturbations we develop a non-linear formalism that 
extends the ordinary Quasi-Linear Theory (QLT) that is routinely used for the purpose. This allows 
us to avoid the usual problem of 90 degree scattering and enable our computation of the mean free 
path of cosmic rays. We apply the formalism to the cosmic ray propagation in the galactic halo and 
in the Warm Ionized medium (WIM). In addition, we address the issue of the transport of cosmic 
rays perpendicular to the mean magnetic field and show that the issue of cosmic ray subdiffusion 
(i.e., propagation with retracing the trajectories backwards, which slows down the diffusion) is only 
important for restricted cases when the ambient turbulence is far from what numerical simulations 
suggest to us. As a result, this work provides formalism that can be applied for calculating cosmic 
ray propagation in a wide variety of circumstances. 

Subject headings: acceleration of particles-cosmic rays-ISM: magnetic fields-MHD-scattering- 
turbulence 



1. INTRODUCTION 

The propagation and acceleration of cosmic rays (CRs) 
is governed by their interactions with magnetic fields. 
Astrophysical magnetic fields are turbulent and, there- 
fore, the resonant and non-resonant (e.g., transient time 
damping, or TTD) interaction of cosmic rays with MHD 
turbulence is the accepted principal mechanism to scatter 
and isotropize cosmic rays (see Schlickeiser 2002). In ad- 
dition, efficient scattering is essential for the acceleration 
of cosmic rays. For instance, scattering of cosmic rays 
back into the shock is a vital component of the first order 
Fermi acceleration (see Longair 1997). At the same time, 
stochastic acceleration by turbulence is entirely based on 
scattering. 

It is generally accepted that properties of turbulence 
are vital for the correct description of CR propagation. 
Historically, the most widely used model is the model 
composed of slab perturbations and 2D MHD perturba- 
tions (see Bieber, Smith, & Matthaeus 1988). The ad- 
vantage of this empirical model is its simplicity and the 
ability to account for the propagation of CRs in magne- 
tosphere given a proper partition of the energy between 
the two types of modes. 

Numerical simulations (see Cho & Vishniac 2000, 
Maron & Goldreich 2001, Miiller & Biskamp 2000, Cho, 
Lazarian & Vishniac 2002, Cho & Lazarian 2002, 2003, 
see also book of Biskamp 2003, as well as, Cho, Lazar- 
ian & Vishniac (2003) and Elmegreen & Scalo 2004 
for reviews), however, do not show 2D modes, but in- 
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stead show Alfvenic modes that exhibit scale-dependent 
anisotropy consistent with predictions in Goldreich & 
Sridhar (1995, henceforth GS95). The approach in the 
latter work makes productive use of the earlier advances 
in understanding of MHD turbulence, that can be traced 
back to Iroshnikov (1963) and Kraichnan (1965) work 
and the classical work that followed (see Montgometry 
& Turner 1981, Higdon 1984, Montgomery, Brown & 
Matthaeus 1987). 

A careful analysis shows that there is no big gap be- 
tween the Reduced MHD and the GS95 model. In fact, 
it was shown in Lazarian & Vishniac (1999) that the nu- 
merical results in Matthaeus et al. (1998) are consistent 
with GS95 predictions. While particular aspects of the 
GS95 model, e.g., the particular value of the spectral 
index, are the subject of controversies'^ (see Miiller & 
Biskamp 2000, Boldyrev 2005, 2006, Beresnyak & Lazar- 
ian 2006, Gogoberidze 2006, Mason et al. 2007), we think 
that, at present, GS95 model provides a good starting 
for developing models of CR scattering, as was done in 
Chandran (2000), Yan & Lazarian (2002, 2004, hence- 
forth YL02, Paper I, respectively), Brunetti & Lazarian 

To address quantitatively these controversies, we need much 
better numerical resolution. For instance, hydrodynamic turbu- 
lence simulations in Kritsuk et al. (2007) showed that only start- 
ing with the 1028'' numerical cubes the bottleneck effects stop 
dominating the measured spectral slope. While the simulations 
in Kowal & Lazarian (2007) show that the bottleneck is less im- 
portant for their MHD code, the exact value of the spectral slope 
is still uncertain. At the same time, the particular theoretically- 
predicted features of turbulence, for instance, the existence of the 
scale-dependent anisotropy, can be reliably established, while the 

2/3 

exact scaling of this dependence, e.g., like fc|| ~ as in GS95 

model or fe|| ~ kj*, where a* < 2/3 (see Beresnyak &; Lazarian 
2006), also require higher resolution simulations. 
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(2007) etc. In particular, the latter three papers used 
the decomposition of MHD turbulence over Alfven, slow 
and fast modes as in Cho & Lazarian (2003) and identify 
the fast modes as the major source of CR scattering in 
interstellar and intracluster medium. 

However, while the turbulence injected on large scales 
may correspond to GS95 model and its extensions to 
compressible medium (Lithwick & Goldreich 2001, Cho 
& Lazarian 2002, 2003), one should not disregard the 
possibilities of generation of additional perturbations by 
CR themselves. Indeed, the slab Alfvenic perturba- 
tion can be created, e.g., via streaming instability (see 
Wentzel 1974, Cesarsky 1980) or kinetic gyroresonance 
instability (see its application for CR transport in Lazar- 
ian & Beresnyak 2006). These perturbations, that are 
present for a range of CRs energies (e.g., < lOOGeV 
for the instabilities above in ISM) owing to non-linear 
damping arising from ambient turbulence (YL02, Paper 
I, Farmer & Goldreich 2004, Lazarian & Beresnyak 2006), 
should also be incorporated into the comprehensive mod- 
els of CR propagation and acceleration. 

At present, the propagation of the CRs is an advanced 
theory, which makes use both of analytical studies and 
numerical simulations. However, these advances have 
been done within the turbulence paradigm which is be- 
ing changed by the current research in the field. As we 
discussed above, instead of the empirical 2D+slab model 
of turbulence, numerical simulations suggest anisotropic 
Alfvenic modes (an analog of 2D, but not an exact one, 
as the anisotropy changes with the scale involved) + fast 
modes or/and slab modes. This calls for important revi- 
sions of the CR propagation, which is the subject of the 
current paper. 

The perturbations of turbulent magnetic field are usu- 
ally accounted for by direct numerical scattering simu- 
lations (Giacalone & Jokipii 1999, Mace et al 2000, Qin 
at al. 2002) or by quasi-linear theory, QLT (see Jokipii 
1966, Schlickeiser 2002). The problem with direct nu- 
merical simulations of scattering is that the present-day 
MHD simulations have rather limited inertial range. At 
the same time, creating synthetic turbulence data which 
would correspond to scale-dependent anisotropy in re- 
spect to the local magnetic field (which corresponds, e.g., 
to GS95 model) is challenging and has not been practi- 
cally realized, as far as we know. 

While QLT allows easily to treat the CR dynamics in a 
local magnetic field system of reference, a key assumption 
in QLT, that the particle's orbit is unperturbed, makes 
one wonder about the limitations of the approximation. 
Indeed, while QLT provides simple physical insights into 
scattering, it is known to have problems. For instance, 
it fails in treating 90° scattering (see Jones, Birming- 
ham & Kaiser 1973, 1978; V51k 1973, 1975; Owens 1974; 
Goldstein 1976; Felice & Kulsrud 2001) and perpendicu- 
lar transport (see Kota & Jokipii 2000, Matthaeus et al. 
2003). 

Indeed, many attempts have been made to improve 
the QLT and various non-linear theories have been at- 
tempted (sec Dupree 1966, Volk 1973, 1975, Jones, 
Kaiser & Birmingham 1973, Goldstein 1976). Currently 
we observe a surge of interest in finding way to go beyond 
QLT. Those include recently developed nonlinear guiding 
center theory (see Matthaeus et al. 2003), weakly nonlin- 
ear theory (Shalchi et al. 2004), second-order quasilinear 



theory (Shalchi 2005a) (see also Shalchi 2006, Webb et al. 
2006, Qin 2007, Le Roux & Webb 2007). At the same 
time, most of the analysis so far has been confined to 
traditional 2D-|-slab models of MHD turbulence. Follow- 
ing the reasoning above, we think that it is important to 
extend the work to the non-linear treatment of CR scat- 
tering to models MHD turbulence that are supported by 
numerical simulations. 

Propagation of CRs perpendicular to the mean mag- 
netic field is another important problem in which QLT 
encounters serious difficulties. Compound diffusion, re- 
sulting from the convolution of diffusion along the mag- 
netic field line and diffusion of field line perpendicular to 
mean field direction, has been invoked to discuss trans- 
port of cosmic rays in the Milky Way (Getmantsev 1963; 
Lingenfelter, Ramary & Fisk 1971; Allan 1972). The 
role of compound diffusion in the acceleration of CRs at 
quasi-perpendicular shocks were investigated by Duffy et 
al. (1995) and Kirk et al. (1996). 

Indeed, the idea of CR transport in the direction per- 
pendicular to the mean magnetic field being dominated 
by the field line random walk (FLRW, Jokipii 1966, 
Jokipii & Parker 1969, Forman et al. 1974) can be easily 
justified only in a restricted situation where the turbu- 
lence perturbations are small and CRs do not scatter 
backwards to retrace their trajectories. If the latter is 
not true, the particle motions arc subdiffusive, i.e., the 
squared distance diffused growing as not as t but as t", 
a < 1, e.g., Q: = 1/2 (Kota & Jokipu 2000, Mace et 
al 2000, Qin at al. 2002, Shalchi 2005b). If true, this 
could indicate a substantial shift in the paradigm of CR 
transport, a shift that surely dwarfs a modification of 
magnetic turbulence model from the 2D+slab to a more 
simulation-motivated model that we deal here. 

It was also proposed that with substantial transverse 
structure, i.e., transverse displacement of field lines, per- 
pendicular diffusion is recovered (Qin et al 2002). Is it 
the case of the MHD turbulence models we deal with? 

How realistic is the subdiffusion in the presence of tur- 
bulence? The answer for this question apparently de- 
pends on the models of turbulence chosen. In this paper 
we again seek the answer for this question within domain 
of numerically tested models of MHD turbulence. 

There arc three major thrusts of the paper: 

I. Extend QLT by taking into account magnetic mirror- 
ing effect on large scales. 

II. Describe CR propagation in Milky Way (e.g., calcu- 
late CR mean free path for different phases of ISM). 
HI. Address the problem of perpendicular transport of 
CR. 

In what follows, we discuss the cosmic ray transport 
in incompressible turbulence in §2. We shall describe 
the §2.1 dispersion of guiding center of CRs and intro- 
duce the broadened resonance function to replace the S 
function in QLT, following which we shall discuss the 
scattering in strong and weak incompressible turbulence 
respectively in §2.2 and §2.3. Then we shall consider the 
scattering by fast modes in §3 and apply the analysis to 
ISM and get mean free path for different phases of ISM 
(§4). In §5, we shall study the perpendicular transport of 
cosmic rays on both large and small scales. We shall also 
discuss the applicability of the subdiffusion. Discussion 
and summary are provided in §6 and §7 respectively. 
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2. CR TRANSPORT IN INCOMPRESSIBLE TURBULENCE 

2.1. General formalism 

It was demonstrated that scattering by the Alfvenic 
turbulence is substantially suppressed due to its 
anisotropy (Chandran 2000, YL02). On the other hand, 
resonant mirror interaction (so-called transit time damp- 
ing or TTD) can arise from the slow modes (also known 
as the pseudo Alfven modes in incompressible limit) and 
it is not subjected to the suppression from anisotropy. 
One may speculate that TTD is the alternative in this 
case. 

The only requirement for TTD is that the projected 
particle speed is comparable to phase speed of the mag- 
netic field compression^ (v|| ~ i^/k\\) for particles to have 
enough collisions with the moving magnetic mirrors be- 
fore they leak out of them. In QLT, this means that only 
particles with a specific pitch angle can be scattered. For 
the rest of the pitch angles, the interaction is still dom- 
inated by gyroresonance, which efficiency is negligibly 
small for the Alfvenic anisotropic turbulence. The quan- 
titative treatment of this interaction is in YL02, which 
used the empirical magnetic field fluctuation tensor from 
Cho et al. (2002), provided the reduction factor of twenty 
orders of magnitude for CR of 100 GeV energies. This 
leads to very large mean free path of CRs. With the res- 
onance broadening, however, we expect that wider range 
of pitch angle can be scattered through TTD, including 
90°. 

The basic assumption of the quasi-linear theory is that 
particles follow unperturbed orbits. In reality, particle's 
pitch angle varies gradually with the variation of the 
magnetic field due to conservation of adiabatic invari- 
ant v\/B, where B is the total strength of the magnetic 
field (see Landau & Lifshits 1975). Since B is varying in 
turbulent field, so is the projection of the particle speed 
and v\\. This results in broadening of the resonance. 
Indeed, the average uncertainty of parallel speed At)|j is 
given by (see Volk 1975) 
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The variation of the velocity is mainly caused by the mag- 
netic perturbation SB^^ in the parallel direction. This is 
true even for the incompressible turbulence we discus- 
sion in this section. For the incompressible turbulence, 
the parallel perturbation arises from the pseudo- Alfven 
modes. The perpendicular perturbation 6Bj_ is higher 
order effect, which we shall neglect in this paper. 

The propagation of a CR can be described as a combi- 
nation of a motion of its guiding center and CR's motion 
about its guiding center. Because of the dispersion of 
the pitch angle A/i and therefore of the parallel speed 
Aw|| , the guiding center is perturbed about the mean po- 
sition < z >= vfit as they move along the field lines. 
As a result, the perturbation 6B{x, t) that the CRs view 
when moving along the field gets a different time depen- 
dence. The characteristic phase function e*'^"^'-*-' of the 
perturbation 6B{x, t) deviates from that for plane waves. 



Assuming the guiding center has a Gaussian distribution 
along the field line, one obtains 



Integrating over z, one gets 
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From Eq.([T|), we obtain 



al Aujj >t^ 



< SBj >\ 



(2) 
(3) 

(4) 



Perturbation 5B\\ exists owing to the pseudo- Alfven 
modes in the incompressible turbulence. Insert the 
Eq.Q into the expression of D^^ (see Volk 1975, Pa- 
per I), we obtain 
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Following are the definitions of the parameters in the 
above equation. Q, fi are the Larmor frequency and pitch 
angle cosine of the CRs. J„ represents Bessel function, 
and w = k±v±/fl ~ k±LRy^l — y?, where R — v/{fll) 
is the dimensionless rigidity of the CRs, L is the injec- 
tion scale of the turbulence. fc_L,A;|| are the components 
of the wave vector k perpendicular and parallel to the 
mean magnetic field, uj is the wave frequency, /"^(k) 
is the energy spectrum of the Alfven modes and I^' (k) 
represents the energy spectrum of magnetosonic modes, 
which in our case at hand are the pseudo- Alfven modes. 
In QLT, the resonance function i?„ = 7r5(fc||W|| —udznfl). 
Now due to the perturbation of the orbit, it should be 

i?„(fc||W|| — UJ ± nfl) 

1/4 

/■oo i{k«v::+nn-uj)t-^kf.v^t^ i 

(1) =^ dte 
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■ exp 



(fc||f /i — UJ -\- nVL)"^ 



exp 



kl^vl 



(6) 



where Ma = SV/va = SB/Bq is the Alfvenic Mach num- 
ber and VA is the Alfven speed. We stress that Eqs. ()5I6|) 
are generic, and applicable to both incompressible and 
compressible medium. 

For gyroresonance [n — ±1,2,...), apparently the re- 
sult is similar to that from QLT for /x ^ A/i = Av\\/v. 
In this limit, Eq.([5]) represents a sharp resonance and be- 
comes equivalent to a (5- function when put into Eq.([5]). 
In general, the result is different from that of QLT, espe- 
cially at a — > 90°, the resonance peak happens at fc||^res ~ 



fl/Av in contrast to the QLT result fcy 



The compressions of the magnetic field in incompressible tur- 
bulence are related the pseudo-Alfven mode, which is the limiting 
case of the slow mode (see Alfven & Falthammar 1963). 



^ We shall show below, that due to the anisotropy, the scattering 
coefficient -D^j^ is still very small if the Alfven and the pseudo- 
Alfven modes are concerned. 
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On the other hand, the dispersion of the W|| means 
that CRs with a much wider range of pitch angle can 
be scattered by the pseudo-Alfven modes through TTD 
(n = 0), which is marginally affected by the anisotropy 
and much more efhcient than the gyroresonance. Below 
we shall consider both the cases for scattering in strong 
and weak turbulence. 

The nonlinear approach we use here is based on parti- 
cle trapping by large scale magnetic perturbations (Volk 
1973,1975). The difference is that we have a Gaussian 
profile (EqO resonance and he adopted a Heaviside step 
function. Formally our approach also has a similarity to 
the second order quasilinear theory that Shalchi (2005a) 
proposed for slab modes although his approach is based 
on a different set of approximations. 

2.2. Strong MHD turbulence 

In strong MHD turbulence, we assume that the Alfven 
and the pseudo-Alfven modes follow the scaling obtained 
in Cho et al. (2002), which is consistent with the GS95 
model: 



Scattering in strong incompressible turbulence 



exp 



Ll/3|fc||| ^ 



(7) 



The pitch angle scattering arising from TTD with 
pseudo-Alfven modes is given by Eq([5]) with n = 0: 
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exp 
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(8) 



where x — kL,x\\ = k\\L are the normalized wavenum- 
bers. limax is the maximum wave vector of the turbu- 
lence, corresponding to the dissipation scale fcmaa;- Since 
all scales contribute to TTD and the magnetic perturba- 
tion increases with scale, the interaction is dominated 
by the large scale moving mirrors (fc ~ 1/^)- Ac- 
cordingly we can make an estimate of the above equa- 
tion. For CRs with small rigidities R = v/{LU) <C 1, 

w = x±Ry^l — fi^ < 1, and Ji{w) ~ w/2. On the large 
scale, the anisotropy of the turbulence is small, and we 
approximate {x\ + in the above expression by 2x\. 

Then we get from Eq.® 
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-(1-^2)5 lE.iqxn 



exp 



(/i - va/v) 
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A^2 



(9) 



where Ei represents the function of exponential integral 
-^1(2;) = dtexp{-xt)/t and q = {x±^maxM^)-^/^ . 
As shown in Figll] the above expression provides a good 
approximation for Eq.®. 

For gyroresonance, the interaction is dominated by the 
first order harmonics because higher order harmonics op- 
erate on smaller scales where turbulence energy is de- 
creasingly small according to the power law spectrum. 
The scattering with the Alfven modes provides 
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Fig. 1. — Scattering of CRs with R = r^/L = 10~* in the strong 
incompressible turbulence [Ma — !)• Solid line represents TTD 
and 'x' line is its analytical approximation (Eq[9]l. Dashdot refers 
to the gyroresonance with the Alfvenic turbulence while dash line 
is for the gyroresonance with the slow modes. 
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Fig. 2.— Scattering of CRs with R = rz,/L = IQ-"' by gyrores- 
onance with the strong Alfvenic turbulence. The solid line repre- 
sents the solution from nonlinear theory (Eg llOI I. and the dashed 
line is the result from QLT (see YL02). 
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Eqs.((8lfT0)l can be evaluated numerically. The results for 
the TTD scattering and gyroresonance are displayed in 
Fig lll2l We see that indeed taking into account the non- 
linear perturbations, CRs of a wide range of pitch angles 
(including 90°) can be scattered through TTD interac- 
tion. At small pitch angles, TTD scattering is negligi- 
ble and only gyroresonance operates. The result for gy- 
roresonance is comparable to that from QLT (see Figl2]), 
showing again that gyroresonance is suppressed for the 
fcj_ 3> fc|| turbulence. 

2.3. Weak MHD turbulence 

In weak turbulence, smaller structure is developed only 
in the perpendicular direction. This entails an increase 
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of k± while keeping fc|| = L^^ unchanged. This type of 
cascade proceeds till "critical balance", — k±v±, is 
reached at the scale Itr- The scaling of the weak turbu- 
lence is Ek oc fcj^ (see Lazarian & Vishniac 1999, ap- 

This renders 

(11) 
In other 



Scattering in weak turbulence 



pendix A there, Galtier et al. 2000). 



LMi 



and 



vtr ^ vlMa = vaM\, where Ma < 1. 
words, the weak turbulence only exists for a limited in- 
ertial range < k± < (LM^)~^ (see also discussion 
in Lazarian 2006). 

If similar to the case of strong MHD turbulence, the 
pseudo-Alfven modes in the weak turbulence follow the 
same scaling as the shear Alfven modes, then Eq.® pro- 
vides 
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exp 
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(12) 



In the case the rigidity R is much less than 1, the Bessel 
function Ji{w) can be replaced by the first order ap- 
proximation /A. The above integration can be then 
evaluated analytically, 
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Fig. 3. — Scattering coefficients in weak incompressible turbu- 
lence. Dashed line represents TTD (i.e., n = 0) in the weak tur- 
bulence on large scales. Solid line and 'x' line (analytical approx- 
imation from Ea ll3l l refer to TTD with the weak pseudo Alfven 
modes on large scales and dashed line represents TTD with the 
strong turbulence on small scales. We see that gyroresonance (i.e., 
n not equal to zero) which is denoted by the dashdot line is negli- 
gible because of the strong anisotropy. 
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-3 Scattering by fast modes (R=1 ' 
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(13) 



The result is compared with the numerical evaluation 
of Eq. in FigO FigO presents the scattering cocfh- 
cients owing to various interactions in an incompressible 
medium with Ma = 0.1. We account for TTD inter- 
actions with both the weak turbulence on large scales 
{L > > Itr) and the strong turbulence on small 

scales (fc^^ < Itr)- Here TTD is present for a smaller 
range of compared to the case of Ma — 1 (Fig[T]), as 

1/2 

Afj, ~ M^' <C fi- For the rest of the pitch angles, the 
turbulence can only scatter CRs through gyroresonance, 
which is inefficient because of the turbulence anisotropy. 



3. CR SCATTERING IN COMPRESSIBLE MHD 
TURBULENCE 

In compressible turbulence, the pseudo-Alfven modes 
become the slow modes, and additional type of pertur- 
bation, fast modes are present. The latter were iden- 
tified as the major scattering agent for the MHD tur- 
bulence that is injected at large scales (YL02,04). The 
papers were using QLT, however. Therefore it is nec- 
essary to check the validity of this conclusion using the 
modified resonance function given by Eq.(IS]). Therefore 
below we provide calculations for scattering induced by 
the fast modes. The scatterings by the Alfven and the 
slow modes are similar to those by the Alfven and the 
pseudo-Alfven modes discussed in the previous section. 
The latter is the similarities of the former claimed on the- 
oretical grounds (GS95) and confirmed using numerical 
simulations (CL03). 

For gyroresonance, according to Ea. (|5l6p . the pitch an- 
gle diffusion coefficient is given by 
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5 Scattering by fast modes (R=1 
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Fig. 4. — Scattering of CRs in the fast mode t urb ulence. The 
solid line is the result for TTD based on NLT (Eq[T5ll and the 'x' 
line represents its analytical approximation (Eq |16l l: the dashdot 
(NLT, Eq[T4ll and dashed (QLT, Paper I) lines refer to the results 
of gyroresonance. Upper: scattering coefficient for CRs of i? = 
rz^/L = 10"*; Lower, scattering coefficient for CRs of J? = r^/L = 
10~^ . This rigidity is enough small that Eq. lll6l l provides a good 
approximation. 
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where ^ = cos 6 is the cosine of wave pitch angle 9. 
For TTD, 
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For CRs with sufhcient smah rigidities R < 1/ (kmaxL), 
the Bessel function can be approximated by the first or- 
der asymptotics and we obtain: 



-v{l — 2 exp 



- VA/vf 



A^2 



(16) 



where kmaxiC) is the cut-off wave number at ^. 

The scattering by the fast modes is exhibited in 
Fig[31 where we provided the numerical evaluations of 
Eqs p4ll5p as well as the analytical appoximation given 
by Ea. p6|) . We adopt viscous damping for the illustra- 
tive calculations (with the physical parameters in WIM, 
see table 1). Realistic calculation for interstellar medium 
will be given in the next section. As we see, the TTD 
interaction dominates for large pitch angles till 90°. Gy- 
roresonance is important for small pitch angles. While 
the TTD interaction is expanded to a much wider range 
(including 90°) compared to QLT result (see, e.g, Fig. 3 in 
Paper I) , the gyroresonance in QLT and nonlinear theory 
(NLT) are comparable. This is because gyroresonance 
happens on a local scale k'^^'^ - unlike TTD. The 
influence of large scale trapping is thus limited. Espe- 
cially for the range ^ > A^, we see marginal difference. 
For large pitch angles, indeed there is a discrepancy be- 
tween NLT and QLT. In fact, we see similar trend in 
the result of Shalchi (2005a), which has close relations 
with the nonlinear theory by Volk (1973, 1975). How- 
ever, because of the dominance of TTD in this range, 
this difference does not count. All in all, nonlinear ef- 
fect is important for CR scattering, particularly for the 
contribution from TTD interaction; for gyroresonance, 
nevertheless, one can use the quasi-linear approximation 
to calculate the scattering of CRs at small pitch angles. 

Note that although the overall contribution from gy- 
roresonance is smaller than that from TTD, gyroreso- 
nance plays an important role in confining the CRs at 
small pitch angles. Without sufficient scattering by gy- 
roresonance (e.g., with the incompressible turbulence, see 
Fig[T]), the mean free path would be unrealistically large 
as TTD is inefficient for the scattering of CRs propagat- 
ing at small pitch angles. Below in the next section we 
shall study the confinement of CRs in the Galaxy by the 
fast modes. 

4. COSMIC RAY PROPAGATION IN GALAXY 

The scattering by fast modes is influenced by the 
medium properties as the fast modes are subject to linear 
damping, e.g.. Landau damping (see Ginzburg 1961). In 
Paper I we showed that the CR scattering is different in 
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Fig. 5. — The turbulence truncation scales in Galactic halo and 
warm ionized medium (WIM). The damping curves flattens around 
90° due to field line wandering (dotted lines, see Paper I, Lazarian, 
Vishniac & Cho 2004); For WIM, both viscous and collisionless 
damping are applicable. 
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Fig. 6. — Scattering efficiency is reduced for moderate energy 
(< ITeV) CRs due to the average of the damping scale around 
90° caused by the field line wandering. The thick lines refer to the 
results without field line wandering while the thin lines represent 
the results taking into account field line wandering. 

different ISM phases, but could not calculate the mean 
free path as we faced the scattering at 90° problem. Us- 
ing the approach above we revisit the problem of the CR 
propagation in the selected phases of the ISM. In par- 
ticular, we shall make quantitative predictions for the 
parallel mean free path of CRs in the galactic Halo and 
Warm Ionized medium (see Table 1 for a list of fiducial 
parameters appropriate for the idealized phases^) assum- 
ing that the turbulence is injected at large scales. 

4.1. Halo 

In Galactic halo (see Table 1), the Coulomb coUisional 
mean free path is ^ lOpc, the plasma is thus in a colli- 
sionless regime. The cascading rate of the fast modes is 

^ The parameters of idealised interstellar phases are a subject of 
debate. Recently, even the entire concept of the phase being stable 
entities has been challenged (see Gazol et al. 2007 and ref. therein). 
Similarly as in Paper I, we were guided in choosing the numbers by 
our communications with Don Cox (2006, private communication). 
However, we accept that different parts of interstellar medium can 
exhibit variations of these parameters (see Wolfire et al. 2003 and 
ref. therein). 
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ISM 


T(K) 


C5(km/s) 


n{cm 


'm/p(cm) 


L(pc) 


B(mG) 


/3 


damping 


Halo 




91 


10"=* 


4 X 1019 


100 


5 


0.14 


Collisionless 


WIM 


8000 


8.1 


0.1 


6 X 10^^ 


50 


6 


0.077 


CoUisionless & viscous 



TABLE 1 

The fiducial parameters of idealized ISM phases and the relevant dampings. The dominant damping mechanism for the fast 

MODES TURBULENCE IS GIVEN IN THE LAST LINE. WIM=WARM IONIZED MEDIUM. 



(Cho & Lazarian 2002). 

r^' = {k/Lf'HV^/V,u (17) 
By equating it with the coUisionless damping rate 



Pitch angle diffusion of CRs in halo 




exp 



(3mi cos^ 6 



(18) 



we obtain the turbulence truncation scale kc = kmax- 

2me 



krL 



4:M\mi cos^ i 
irrrieP sin** 9 



exp 



firrii cos^ 9 



(19) 



where [3 = Pgas/Pmag- 

The scale kc depends on the wave pitch angle 9, which 
makes the damping anisotropic. As the turbulence un- 
dergoes turbulent cascade or/and the waves propagate 
in a turbulent medium, the angle 9 is changing. As dis- 
cussed in Paper I the field wandering defines the spread 
of angles. During one cascading time, the fast modes 
propagate a distance VTcas and see an angular deviation 

tsai69 : 



tan^ 



'II 



tan^ 



which is 



tan; 



M\ cos 9 




27(fcL)i/2 + 


V kL ) 



(20) 



As evident, the damping scale given by Eq. (|T9)) varies 
considerably especially when — > and 9 — + 90°. For the 
quasi-parallel modes, the randomization (oc (fcL)~^/^) is 
negligible since the turbulence cascade continues to very 
small scales. On small scales, most energy of the fast 
modes is contained in these quasi-parallel modes (Paper 
I, Petrosian, Yan & Lazrian 2006). 

For the quasi-perpendicular modes, the damping rate 
(Eq fT^ should be averaged over the range 90° — 50 ~ 90°. 
Equating Eq. fTTl) and Eq. (fT8|) averaged over 59, we get 
the averaged damping wave number (see FiglS]). The 
field line wandering has a marginal effect on the gyrores- 
onance, whose interaction with the quasi-perpendicular 
modes is negligible (Paper I). However, TTD scattering 
rates of moderate energy CRs (< lOTeV) will be de- 
creased owing to the increase of the damping around the 
90° (see Fig|6]). For higher energy CRs, the influence of 
damping is marginal and so is that of field line wander- 
ing. 

We adopt QLT for calculating the gyroresonance (see 
Paper I). As we see from FigHl the QLT result in the 
range /i > A/i provides a good approximation for our 
calculations using the non-linear approximation given by 
Eq. (fT4|) . For CRs with sufficiently small rigidities, the 
resonant fast modes (fc^es ~ l/(^/^)) are on small scales 
with a quasi-slab structure (see Fig[S]). For the scattering 
by these quasi-parallel modes, the analytical result that 
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Fig. 7. — Pitch angle diffusion coefficients in halo. Upper lines 
in the plots represent the contribution from TTD and lower lines 
are for gyroresonance. 
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Pitch angle diffusion of CRs in WIM 
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Fig. 8. — Pitch angle diffusion coefficients in WIM. The notations 
are the same as in Fig[7] 

follows from QLT approximation (see Paper I) for the 
gyroresonance is'' 



D 



pp 



7rt;/i°-5(l-^2^ 



21 



\\l + {R^ifri - (tan2 0,-Kl)-5 
^ {[1 + {Riifr 'i - (tan2 9, + 1)"! 

where tan 0c = fc_L_c/fc||,res- This justifies our use of the 
analytical approximation above. 

Once we the know the functional form of the D^^, we 
can obtain the corresponding mean free path (Earl 1974): 



4./0 



(22) 



^ It can be shown that the QLT result follows from our more 
general result Ea. lll4l l if we put A/i ^ 0. 
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Confinement of CRs in different pfiases of ISM 




where Dj,,, canbe obtained by Eqs. (|15IIIBl [TO]) . Insert 
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Fig. 9. — The mean free paths in two different phases of ISM: 
halo (solid hne) and WIM (dashed hne). At lower energies (< 
lOOGeV), the different dependence in WIM is owing to the viscous 
damping (sge Fig[5]and text J. 

DJ^^ and QLT result for D'^^ into the above expression, 
we get the mean free path of CRs in halo. 

The mean free path is sensitive to the scattering by gy- 
roresonance at small pitch angles, due to the influence of 
damping on the fast modes on small scales. Fig [7| shows 
the pitch angle diffusion of CRs with different energies 
due to the TTD and gyroresonance. 

The weak dependence of the mean free path (see FigE]) 
of the moderate energy (e.g< ITeV) CRs in halo results 
from the fact that gyroresonance changes marginally 
with the CR energy (see Figl7|). Gyroresonance happens 
on small scales where the fast modes develop a quasi-slab 
structure because of the damping (see Fig[5]). In the case 
of halo, the critical 9c changes more slowly compared to 
the case in WIM (see FiglH]), the scattering by gyrores- 
onance is thus marginally changing with the energy (see 
Eql2T] and Figl7|). For higher energy CRs with larger gy- 
roscales, the influence of damping is small, and thus the 
CR mean free path begins increasing with energy. 

4.2. Warm Ionized Medium 

In warm ionized medium, the Coulomb coUisional 
mean free path is Imfp = 6 x lO^^cm and the plasma 
P ~ 0.11. Suppose that the turbulence energy is in- 
jected from large scale, then the compressible turbulence 
is subjected to the viscous damping besides the colli- 
sionless damping. By equating the viscous damping rate 
with the cascading rate (Eq[T7|) , we obtain the following 
truncation scale. 



where Xc — ^I^^J" j '7o is the longitudinal viscosity. 
In the low (3 regime, the motions are primarily perpendic- 
ular to the magnetic field so that dvx/dx = n/n ~ B/B. 
The longitudinal viscosity enters here as the result of 
distortion of the Maxwellian distribution (see Bragin- 
skii 1965). The transverse energy of the ions increases 
during compression because of the conservation of adia- 
batic invariant v'^/B. If the rate of compression is faster 
than that of collisions, the ion distribution in the momen- 
tum space is bound to be distorted from the Maxwellian 



(i-e2)-i/3«i 



(23) 



isotropic sphere to an oblate spheroid with the long axis 
perpendicular to the magnetic field. As a result, the 
transverse pressure gets greater than the longitudinal 
pressure, resulting in a stress ~ rjodv^/dx. The restora- 
tion of the equilibrium increases the entropy and causes 
the dissipation of energy. 

The viscous damping scale is compared to coUisionless 
cutoff scale fEalT9| in Fig[5] As shown there, both vis- 
cous damping and coUisionless damping are important 
in WIM. Viscous damping is dominant for small 6 and 
coUisionless damping takes over for large 9 except for 
9 = 90°. This is because coUisionless damping increases 
with 9 much faster than the viscous damping. For suf- 
ficiently small wave pitch angles, the viscous damping is 
too small to prevent the fast modes to cascade down to 
scales smaller than the mean free path Imfp- Because of 
the similar quasi-slab structure on small scales, Ea. (PT|) 
can be also applied in WIM. The results are illustrated 
in Fig[5] Compared to the case in halo, we see that the 
qualitative difference stands in the gyroresonance. This 
is because gyroresonance is sensitive to the quasi-slab 
modes for which dampings differ in halo and WIM. 

The mean free paths of CRs are given in Figl9l We 
see the mean free path first decreases with the energy 
till lOOGeV. This is because of the influence of viscous 
damping on gyroresonance (see Fig [5]). The lower the 
energy of the CRs is, the narrower is the wave vector cone 
of the available fast modes and thus the less efficient is 
the gyroresonance according to Eq. ([2T|) (see FigE]). For 
higher energy CRs, of which rigidities are larger than R > 
(3 X 10^)^^, the maximum wave number of the resonant 
wave modes are determined by coUisionless damping (see 
Fig[5|). As a result, the mean free path grows with energy 
similar to the case in halo where coUisionless damping is 
dominant. 

4.3. Other phases 

In hot ionized medium (HIM), the plasma is also in 
coUisionless regime, but the density is higher and the 
plasma beta is larger than 1. The damping by protons 
thus becomes substantial especially at small pitch angles. 
The damping truncates the turbulence at much larger 
scales than the gyroscales of the CRs of the energy range 
we consider. No gyroresonance can happen and some 
other mechanisms are necessary to prevent CRs stream- 
ing freely along the field. The turbulence injected from 
small scales might play an important role (see §5). 

In partially ionized gas one should take into account an 
additional damping that arises from ion-neutral collisions 
(see Kulsrud & Pearce 1969, Lithwick & Goldreich 2001, 
Lazarian, Vishniac & Cho 2004). In the latter work a 
viscosity-damped regime of turbulence was predicted at 
scales less the scale k'^l^j^f, at which the ordinary mag- 
netic turbulence is damped by ionic viscosity. The corre- 
sponding numerical work, e.g., Cho, Lazarian & Vishniac 
(2003) testifies that for the viscosity-damped regime the 
parallel scale stays equal to the scale of the ambipolar 



damping, i.e., ki\ 



(,, while k± increases. In that 



respect, the scattering by such magnetic fluctuations is 
analogous to the scattering induced by the weak turbu- 
lence (see §2.3). The difference stems from the spectrum 
of k± is shallower than the spectrum of the weak tur- 
bulence. The predicted values of the spectrum for the 
viscosity-damped turbulence E{k±^) ^ kj^^ (Lazarian et 
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al. 2004) are in rough agreement with simulations. More 
detailed studies of scattering in partially ionized gas will 
be provided elsewhere. 

5. PERPENDICULAR TRANSPORT 

While in the earlier sections we dealt entirely with the 
diffusion parallel to the magnetic field, in this section we 
deal with the diffusion perpendicular to the mean mag- 
netic field. The assumption that CRs follow the magnetic 
field averaged over their Larmor radius is pretty accurate 
in most situations, e.g., for Galactic CRs. 

Compound diffusion happens when particles are re- 
stricted to the magnetic field lines and perpendicular 
transport is solely due to the random walk of field line 
wandering (see Kota & Jokipii 2000). In the three- 
dimensional turbulence, field lines are diverging away 
due to shearing by the Alfven modes (see Lazarian & 
Vishniac 1999, Narayan & Medvedev 2002, Lazarian 
2006, 2007). Since the Larmor radii of CRs are much 
larger than the minimum scale of eddies l±^min, field 
lines within the CR Larmor orbit are effectively diverg- 
ing away owing to shear by the Alfvenic turbulence. The 
cross-field transport thus results from the deviations of 
field lines at small scales, as well as field line random 
walk at large scale (> min[_L/M^, L]). 

Most recently the diffusion in magnetic fields was con- 
sidered for thermal particles in Lazarian (2006, 2007). 
In what follows we modify the results of these studies for 
the case of CRs. 

For perpendicular diffusion, the important issue is the 
reference of frame. We emphasize that we consider the 
diffusion perpendicular to the mean field direction in the 
global reference of frame. 

5.1. Perpendicular diffusion on large scale 

High Ma turbulence: 
High Ma turbulence corresponds to the field that is eas- 
ily bended by the hydrodynamic motions at the injection 
scale as the hydro energy at the injection scale is much 
larger than the magnetic energy, i.e. pV^ ^ . The 
turbulence in clusters of galaxies is the high Ma turbu- 
lence. In this case the magnetic field becomes dynami- 
cally important on a much smaller scale, i.e., the scale 
Ia = L/M\ (see Lazarian 2006). If A|| » Ia, the CR dif- 
fusion is controlled by the straightness of the field lines, 
and 

D^^D\\K 1/^Iav, Ma > 1, A|| > Ia. (24) 

The diffusion is isotropic if scales larger than I a are con- 
cerned. 

In the opposite limit A|| < Ia, the stiffness of B field is 
negligible. The CR diffusion is insensitive to the topology 
of the magnetic field. In the global reference of frame, 
there is no distinction between the perpendicular and 
parallel direction. Naturally, a result for isotropic turbu- 
lence, namely, 



Di =Di 



l/3A||t;, 



(25) 



holds. 

Low Ma turbulence: 
For strong magnetic field, i.e., the field that cannot be 
easily bended at the turbulence injection scale, individual 
magnetic field lines are aligned with the mean magnetic 
field. The diffusion in this case is anisotropic. As we 



mentioned earlier, if the turbulence is injected at scale 
L it stays weak for the scales larger than Itr given by 
Eq. pT|) and it is strong at smaller scales. 

Consider first the case of the CR parallel mean free 
path larger than the injection scale of the turbulence, i.e., 
A|| > L. The perturbations of the field are uncorrelated 
over scales larger than LM'^ in the direction perpendic- 
ular to the mean magnetic field. Indeed, this perpendic- 
ular distance corresponds to the particle moving parallel 
distance of the order L, which is the scale of the energy 
injection, which in a simplified picture of turbulence^ is 
the maximal scale over which the magnetic perturbations 
are correlated. 

In this situation the random walk steps in perpendic- 
ular direction are of Itr length. Thus, to diffuse over a 
distance R with random walk of Itr one requires {R/ltr)"^ 
steps. The time of the individual step is L/v\\, then 



6t {R/krYL/vw 



l/iLvM\, Ma < 1, A|| > L. 

(26) 



This is similar to the case discussed in the FLRW model 
(Jokipii 1966). However, we obtain the dependence of 
M\ instead of their M\ scaling. This difference is not 
crucial for the environment like Milky Way or solar wind, 
for which Ma ~ 1, but may be important for other envi- 
ronments where strong slightly perturbed magnetic field 
is present, e.g., solar corona. 

What would be the CR diffusion perpendicular to the 
mean magnetic field in the opposite case of A|| < LI 
The time of the individual step is L'^ /D\\. Therefore the 
perpendicular diffusion coefficient is 

i?2 i?2 



D 



6t 



{Rlltr?L^lD\\ 



Ma < 1, A|| < L 



(27) 

which coincides with the result obtained for the diffusion 
of electrons in magnetized plasma (Lazarian 2006). The 
turbulence in the interplanetary medium is in this regime 
with AIa < 1. From the Eq. ([771) . we obtain a constant 
ratio of A^/A|| = D±/D^^ — M\, consistent with the 
Palmer consensus (Palmer 1982). 

We mention parenthetically, that our arguments above 
can be repeated for any random walk process in the per- 
pendicular direction with a step Sx. If we can write 
D± = {Sx/Sz^D^^, we, naturally, recover the result in 
Eq. (HIl). 

5.2. Perpendicular diffusion on small scales 

The diffusion of CR on the scales <C L may be different. 
We consider particular examples below. 

High Ma turbulence: 
Consider the diffusion on scales that are On scales Ay < 

fc^^ < Ia, i.e., on scales at which CR are in diffusive 
regime, but the magnetic fields are strong enough to in- 
fiuence turbulent motions, the mean deviation of a field 
in a distance fcj^^ = 6z is given by (Lazarian & Vishniac 

1999), 



< {Sxf >i/2= 



([fe]MA)3/2 
33/2^1/2 ■' 



Ma > 1 (28) 



° For the sake of simplicity we disregard the effects of the in- 
verse cascade that can increase the correlation scale of magnetic 
perturbations. 
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Thus, for scales much less than L 

(29) 

which for a limiting case fc|| ^ 1^ gets the result consis- 
tent with Eq. ((M)) . 

Low Ma turbulence: 
On scales larger than Itr, the turbulence is weak (see 
§3.2). The mean deviation of a field in a distance 6z is 
given by Lazarian (2006): 

< (fe)2 >i/2^ _|f]!^!_Ml, Ma<1. (30) 

For the scales L > k\^^ = 5z > A|| we combine Eq. (|30p 
with 

5z = y^^ii^t (31) 
and get for scales much less than L 

" ^ = ^^'^ ^ (32) 
which for a limiting case of ^ L^^ coincides up to a 
factor with the Eq. (gT]). 

Eqs. (p9| and ((32| certify that the perpendicular dif- 
fusion at scales much less than the injection scale accel- 
erates as z grows. The reason is that there is no random 
walk on small scales up to the injection scale of the strong 
MHD turbulence {Itr for Ma < 1 and I a for > 1). 
The diffusion on scales less than the turbulent injection 
scale is important for describing propagation and accel- 
eration of CR in supernovae shells, clusters of galaxies 
etc. 

5.3. Subdiffusion 

The diffusion coefhcient in Eq. (gT]), i.e., £»||M|, means 
that the transport perpendicular to the dynamically 
strong magnetic field is a diffusion, rather than subd- 
iffusion, as it was stated in a number of recent papers. 
Let us clarify this point by obtaining the necessary con- 
ditions for the subdiffusion to take place. 

In the papers discussing compound diffusion (see Kota 
& Jokipii 2000, Webb et al. 2006), iSx)^/Sz = D^pat is a 
spatial diffusion constant coefficient. If we adopt this, we 
shall indeed get from Eq. (|3ip the perpendicular diffusion 
coefficient 

= (^^^ = DspatD\i/Sz = D,patD\'\5t)-^/^ 

(,33) 

Therefore the perpendicular transposition will be x — 

D±St = DspatD^\'^{5tY^'^ in accordance with the findings 

in the aforementioned papers. 

The major implicit assumption in the reasoning above 
is that the particles trace back their trajectories in x 
direction on the scale 5z. If this is not true, the intro- 
duction of the diffusion coefficient Dgpat does not make 
sense. 

When is it possible to talk about tracing particle tra- 
jectories back? In the case of random motions at a single 
scale only, the distance over which the particle trajecto- 
ries get uncorrelated is given by the Rechester & Rosen- 
bluth (1978) model. On scales larger than the Rechester 



& Rosenbluth scale fc|| ^ > Ljir, the separation between 

field lines 5x grows monochromatically with the distance 
5z along the magnetic field, no retracing can happen in 
this case Assuming that the damping scale of the turbu- 
lence is larger than the CR Larmor radius, this model, 
when generalized to the anisotropic turbulence provides 
(Narayan & Medvedev 2001, Lazarian 2006) 

Lrr = l\\ ln(/±,mm/rL) (34) 

where l\\^rnin is the parallel scale of the cut-off of turbu- 
lent motions, l±^rnin is the corresponding perpendicular 
scale, is the CR Larmor radius. The assumption of 
< l±,min can be valid, for instance, for the Alfvenic 
motions in partially ionized gas. However, it is easy to 
see that, even in this case, the corresponding scale is 
rather small and therefore subdiffusion is not applicable 
for the transport of particles in the Alfvenic turbulence 
over scales ^ l\\,min- 

If r/, > l±^mim as it is a usual case for the Alfvcn mo- 
tions in the phase of ISM with the ionization larger than 
« 93%, where the Alfvenic motions go to the thermal 
particle gyroradius (see the estimates in Lithwick & Gol- 
dreich 2001, Lazarian et al. 2004), the subdiffusion of CR 
is not an applicable concept for the Alfvenic turbulence. 
This does not preclude subdiffusion from taking place 
in particular models of magnetic perturbations, e.g., in 
the slab model considered in Shalchi (2005b), but we be- 
lieve in the omnipresence of the Alfvenic turbulence in 
interstellar gas (see Armstrong et al. 1995). 

6. DISCUSSIONS 

The present paper extends our study in Paper I. As 
in Paper I we mostly deal with the magnetic perturba- 
tions that are part of the large scale turbulent cascade, 
which is consistent with the Big Power Law in the sky 
observed via radio-scattering and scintillation technique 
(Armstrong et al. 1995). In both papers we use the 
description of the MHD turbulence that follows from nu- 
merical simulations. 

In Paper I we have the CR scattering calculated in the 
selected interstellar environments making use of Quasi- 
Linear Theory (QLT). Because of the limitations of the 
QLT, we could not provide calculations of the mean free 
path in Paper I, which limited the utility of the study. 
In this paper we extended the non-linear approach sug- 
gested in Volk (1975) to treat the scattering, which al- 
lows us to calculate the mean free paths that arise from 
CR interactions with the fast modes. In doing so, simi- 
lar to Paper I, we take into account damping of the fast 
modes in the presence of the field wondering induced by 
the Alfvenic modes. 

Our results show that in WIM and halo of our Galaxy, 
confinement of bulk CRs are mostly due to the compress- 
ible modes. We obtain CR mean free paths about a few 
parsec, consistent with what observations indicate. The 
major difference with earlier picture is the dependencies 
of CR transport parameters on the medium properties. 
The dependence appears as a result of damping of the 
fast modes. For low energy CRs (< lOOGeV), if dom- 
inated by viscous damping, the mean free path of CRs 
would decrease with energy; with collisionless damping, 
however, CRs' mean free path stays almost a constant. 
Field line wandering in general increases the damping of 
the fast modes and reduces the scattering efficiency of 
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CRs. For higher energy CRs, the influence of damping is 
limited, and their mean free path increases with energy. 

The dependencies on the turbulence damping and 
therefore the phase properties should have various im- 
plications from ratio of secondary to primary elements, 
diffuse Galactic 7 ray emission, to the CMB synchrotron 
foreground. With precise measurements, the under- 
standing of CMB is now constrained by our understand- 
ing of the foreground. The variation of CR index over the 
Galaxy may paralyze the synchrotron templates. Such 
variations can be addressed on the basis of the more elab- 
orate CR propagation theory. 

The importance of this study goes beyond the interstel- 
lar medimn. For instance, Brunetti & Lazarian (2007), 
treated acceleration of CRs for plasma in clusters of 
galaxies appealing to the fast modes, which is the ap- 
proach to CRs similar to that in Paper I. We believe that 
the non-linear treatment may be useful for such cases as 
well. In addition, stochastic acceleration by the MHD 
turbulence is a promising mechanism for generating high 
energy particles during solar flares (see, e.g., Petrosian & 
Liu 2004, and references therein). An application to the 
acceleration of CRs in solar flares will be given in Yan, 
Lazarian & Petrosian (2007, in preparation). 

In our treatment we attempted to use the scalings that 
(a) are consistent with numerical calculations and (b) 
whose amplitudes we can estimate with a sufficient de- 
gree of precision. Therefore our present study does not 
deal with scattering of CRs by the fast modes on the 
scales I > LM\, Ma < 1, i.e., on the scales where the 
Alfvenic turbulence in the weak regime. It was suggested 
by Chandran (2005) that the weak fast modes at small 
pitch angles tend to steepen due to the coupling with 
the Alfven modes. When the resulting scaling of the fast 
modes becomes clearer, our approach will be applicable 
to them. 

We have not quantitatively dealt in the present paper 
with the case of the slab Alfven modes created by insta- 
bilities^. The CR scattering by the perturbations created 
by those modes may dominate over the gyroresonance 
with the fast modes, especially for CRs of low energies, 
i.e., whose gyroresonance with the fast modes is ineffl- 
cient due to the fast modes damping (see estimates in 
Lazarian & Beresnyak 2006). Progress in quantitative 

® Streaming instability (see Cesarsky 1980) is an example of 
such instability. However, the instability is suppressed by both 
ion-neutral damping (Kulsrud & Pierce 1969) and the ambient 
turbulence (YL02, Farmer &; Goldreich 2004, Paper I, Lazarian 
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description of the non-linear stages of the instabilities 
that can create slab modes should enable comprehensive 
models that include both the fast modes and the slab 
modes. 

In addition, we addressed the issue of perpendicular 
diffusion, the issue that we have not dealt with in Paper 
I. We found, that similar to the case of thermal diffusion 
discussed in Lazarian (2006), the diffusion of CRs de- 
pends on the Alfvenic Mach number Ma- We found that 
the suppression of the perpendicular diffusion compared 
to the parallel one scales as M\ for Ma < 1. Approach- 
ing the issue of subdiffusion, we found that it is negligible 
for CRs in the Alfvenic turbulence. 

7. SUMMARY 

Our result can be briefly summarized as follows 

• Treatment of the scattering in both strong and 
weak MHD turbulence has been generalized to ac- 
count for perturbations of the particle orbits. We 
found that the non-linear treatment is essential for 
calculating mean free paths of CRs. 

• Our calculations of scattering rates performed for 
different modes of MHD turbulence, assuming that 
the turbulence is injected at large scales, confirm 
the dominance of the fast modes for scattering of 
the bulk CRs in WIM and Galactic halo. 

• We obtained the relation between the CR diffusion 
coefficient parallel to the magnetic field and the CR 
diffusion coefficient perpendicular to the magnetic 
field. We show that CR transport perpendicular to 
the magnetic field depends strongly on the Alfven 
Mach number of the turbulence. 
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